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Abstract 


Nuclear reactions in biological systems produce low-energy frag- 
ments of the target nuclei seen as local events of high linear energy 
transfer (LET). A nuclear-reaction formalism is used to evaluate 
the nuclear-induced fields within biosystems and their effects within 
several biological models . On the basis of direct ionization inter- 
action, one anticipates high-energy protons to have a quality factor 
and relative biological effectiveness (RBE) of unity. Target fragmen- 
tation contributions raise the effective quality factor of 10 GeV pro- 
tons to 3.3 in reasonable agreement with RBE values for induced 
micronuclei in bean sprouts. Application of the Katz model indicates 
that the relative increase in RBE with decreasing exposure observed 
in cell survival experiments with 160 MeV protons is related solely 
to target fragmentation events. Target fragment contributions to lens 
opacity give an RBE of 1.4 for 2 GeV protons in agreement with the 
work of Lett and Cox. Predictions are made for the effective RBE 
for Harderian gland tumors induced by high-energy protons. An ex- 
posure model for lifetime cancer risk is derived from NCRP 98 risk 
tables, and protraction effects are examined for proton and helium 
ion exposures. The implications of dose rate enhancement effects on 
space radiation protection are considered. 

Introduction 

As an energetic-charged particle passes through a region of material, it will suffer many 
atomic/molecular interactions to which only small amounts of energy are given to ioniza- 
tion/excitation at each interaction site. Secondary electrons and photons propagate the energy 
from the initial loss site causing a broadening of the particle track (Katz et al. 1971; Chatterjee, 
Maccabee, and Tobias 1973; Kellerer and Chmelevsky 1975; Paretzke 1988). In this way, the 
passing particle can affect a localized volume, even though the path is remote to the localized 
volume itself. Occasionally, the passing energetic particle undergoes a nuclear reaction in which 
a large amount of its kinetic energy is given to the nucleus of the struck atom. Often, several 
nuclear disintegration fragments (nuclear stars) are produced of sufficient energy to form well- 
defined tracks emanating from the interaction site. These fragments may also affect localized 
volumes remote to the initial particle trajectory. 

In the present work, we endeavor to evaluate the field quantities for the nuclear disintegration 
products formed in nuclear reactions in tissue. We consider explicitly the localized energy deposit 
associated with the fragments and their resulting impact on biological response in the context of 
several dose response models. The effect of secondary electrons is studied through the Katz track 
structure model. All other models used depend only on the density of ionization along the track 
length. Several biological end points are considered including cell survival, cell transformation, 
cataract formation, Harderian gland tumors, and cancer risk. Future work will be oriented 
toward track structure models. 



Target Fragment Transport 

An expression for ion fluence in a region bounded by a surface T (Wilson 1977) is given as 

^ (x ' a B ) “ sj(E) Pj (£)^ a £-,) 

X^{x+ [Rj ( E ) - Rj (£')] n, E") (1) 

where = RJ l [p- d+ Rj(E)}, p=Flx, d = n T n , x , and Pj(E ) is the total nuclear survival 
probability. In equation (1), S(E ) denotes the stopping power, E denotes the energy, f l denotes 
the direction of motion, Rj{E) denotes the range, and Aj denotes the atomic mass number of 
ions of type j. The integral over E' is a summation over the collisional source distribution from 
the boundary {E' = Ej) to the point x(£' = E). We approximate equation (1) in a perturbation 
series by taking 

where 

x' = x + [Rj (E) - Rj (E f )] n p' = n' ■ x', d! = ft ' - r nv , and E% = R- k l [p' -d' + R k (£")] 

We specialize to a unidirectional monoencrgetic beam of type M particles at the boundary as 

<f> k (T,fl,E) = <5 (E - Eo) 6 (n - Oo) (3) 

for which equation (1) may be simplified to 

4>J (x, ft, E) = e-”j(p-Q t (El - Eo) 8 (O - ft 0 ) 

+ [ 1 dE‘ A E- exp {-<Tj [Rj (E') - Rj (E)] [~a rn (p - d') r jM (E r , E", JI, fi„)] } (4) 

JE \E) - ' 

where E" = R7 l [R k (E 0 ) - p' + d']<jj is the nuclear macroscopic cross section and rjM is the 
corresponding fragmentation cross section. If we restrict ourselves to a small volume of material 
(aj(p-d) < 0,01 corresponding to a few to several millimeters), then 


4>j (x, E) = 8 (E- f - E 0 ) 6 (ft - n 0 ) 


+ J* dE 's^E) (£/ ’ E " 5 n + r & 


( 5 ) 


where the projectile and target fragment terms are shown separately. The projectile term has 
a contribution only at energies near E 0 , whereas the target term has a contribution only for 
E E 0 * At high energy we have 


Ei » E + -Sj{E) [p - d) 

A j 


( 6 ) 
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in which case 


4>J (X, n, E) « <5 (E - E 0 ) 6(0- Oo) 

+ r £v/ (-£> o, n 0 ) (p-d) 


The absorbed energy density is then 

Dw =?r je J dn Sj (E) (x, n, e) 


= -s'm (fi'o) + iS j (^o) °"m (p - d ) 
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The first term clearly dominates the second (that is, [<jy/(/9 — d) < 0.01 and 
Sj\j(E 0 ) ;$> ^Sj(i? 0 ) rnj^/] ) • The fact that the third term is nonnegligible results from the 
j 

large stopping power of the low-energy fragments represented by the third term compared with 
S m {E 0 ), which also results in all their energy being deposited locally (Wilson 1977). On this 
basis, equation (7) may be reduced to 


4>j (x, n, E) « S jM 6 (E - Eo) 6 (SI- Slo) 

+ dE 'sy(E) r f M E °’ ^ (9) 

Accordingly, the high-energy beam exposure of a small object can be treated by evaluating the 
direct ionization of the primary particles and the transport of low-energy fragments produced 
uniformly throughout the volume. This treatment is represented in equation (9). We now 
consider some applications of target fragment transport. 

Effects in Conventional Risk Assessment 

Biological risks are related to the local energy deposited by the passage of energetic ions. 
The ionization energy loss, the (Sj[f(E 0 ) term of eq. (8)), is on the order of 0.2 Z 2 keV//im for a 
passing relativistic ion of charge Z . Some ions produce nuclear reactions in which 10 to 100 MeV 
(per event) are released (Wilson, Stith, and Stock 1983) locally as secondary nuclear fragments 
(the third term of eq. (8)). The average energy transfer rate is approximately 0.05A 2 / 3 keV/^m, 
where A is the ion atomic weight. Because the quality factor of the fragments is usually taken 
as 20, the risk associated with direct atomic ionization is on the order of the risk associated 
with the nuclear events for incident low-charge ions (Z < 5), whereas the risk is dominated by 
the direct ionization for high charge (Z 5). At a sufficiently low energy, the direct ionization 
always dominates the biological risk independent of the ion charge and mass. In the present 
section, we quantify these various contributions to biological risk using quality factors presently 
in force (ICRP 26 in Anon. 1987) and evaluate the effects of newly proposed quality factors 
(ICRU 40 in Anon. 1986; ICRP 60 in Anon. 1991). 
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We consider a volume of tissue through which a monoenergetic ion fluence 4> z {E p ) of energy E p 
has passed and evaluate the energy absorbed by the media in the passage. Several processes exist 
by which the ion gives up energy to the media: electronic excitation/ionization, nuclear coulomb 
scattering, nuclear elastic scattering, and nuclear reaction. The electronic excitation/ionization is 
contained in the stopping power that is evaluated by methods discussed in relation to equation (8) 
(Wilson et al. 1984). The nuclear coulomb elastic scattering is highly peaked at low-momentum 
transfer, and the energy transfers per event of a few hundred eV or less are typical (Wilson 
et al. 1983). The nuclear elastic scattering energy transfer is on the order of 1 MeV or less and 
can be neglected in comparison with reactive processes. A model for proton-induced reactions 
in tissue constituents has been given elsewhere (Wilson et al. 1988, 1989) and will provide the 
basis for the present evaluation. 

The secondary-particle radiation fields <j>j(E) are given as 

= S~CE) J E (10) 

where Sj(E) is the stopping power and Cj(E) is the particle source energy distribution, which 
is given as 

Cj (E) =poj (E p ) fj (E) 4> z (E p ) (11) 

where p is the nuclear density, aj(Ep) is the fragmentation cross section, and fj(E ) is the 
fragment spectrum as discussed elsewhere (Wilson, Townsend, and Khan 1989). The total 
absorbed energy is approximately 

poo 

D z (E P ) = S z (E p ) 4> z (E p ) +Y, S 3 (E) 4>j (E) dE (12) 

j Jo 

Equation (12) may be written as 

D z (Ep) = s z (Ep) 4>z (Ep) + Ej <7j (Ep) p <j) Z ( E z ) (13) 

where Ej is the average energy associated with each spectral distribution fj{E). Similarly, the 
dose equivalent is 


Hz (E P ) = Qz (E P ) S z (E p ) <p z (E P ) 

+ Y^QFjEj <Tj (E p ) p 4> z (Ep) (14) 

j 

where Qp^ is the spectral-averaged quality factor of the jth secondary particle (Shinn, Wilson, 
and Ngo 1990). The sum over j will include the usual “evaporation” products, including the 
low-energy protons. 

Evaluation of Dose Equivalent 

We now evaluate equation (14) for the conventional quality factor Q(L ) that is dependent on 
linear energy transfer (LET) (ICRP 26 in Anon. 1987; ICRP 60 in Anon. 1991) and the lineal- 
energy-dependent quality factor Q(y ) that was recently proposed (ICRU 40 in Anon. 1986). To 
implement Q(y ), we used appendix B of ICRU 40 in which the lineal energy distributions are 
assumed to be linearly dependent on y at a fixed LET. Some problems of this assumption have 
been discussed by Townsend, Wilson, and Cucinotta (1987), which we circumvent herein by 
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Table 1. Spectral- Averaged Quality Factor for Individual Isotopes 
Produced in 1-GeV Proton-Induced Reactions in 16 0 


h 


o' 

< 

Qfj 

ICRP 26 

ICRU 40 

ICRP 60 

1 

1 

8.69 

2.73 

3.71 

2.65 

1 

2 

10.70 

4.09 

5.87 

4.36 

1 

3 

10.40 

5.20 

7.72 

6.07 

2 

3 

11.20 

12.38 

19.52 

17.86 

2 

4 

12.30 

13.90 

21.88 

20.16 

3 

6 

6.85 

19.25 

23.19 

20.03 

3 

7 

6.16 

19.50 

22.38 

19.20 

4 

9 

4.79 

19.77 

15.14 

15.06 

4 

10 

4.11 

19.73 

14.95 

15.01 

5 

9 

4.79 

19.81 

11.47 

13.02 

5 

10 

4.11 

19.78 

11.54 

13.06 

5 

11 

3.71 

19.75 

11.88 

13.26 

5 

12 

2.74 

19.70 

12.58 

13.70 

6 

11 

3.71 

19.79 

9.96 

12.03 

6 

12 

2.74 

19.75 

10.71 

12.48 

6 

13 

2.05 

19.68 

11.96 

13.28 

6 

14 

1.34 

19.56 

14.11 

14.53 

6 

15 

0.69 

19.17 

18.30 

17.06 

7 

13 

2.05 

19.74 

10.53 

12.32 

7 

14 

1.37 

19.64 

12.61 

13.62 

7 

15 

0.69 

19.35 

16.73 

16.14 

8 

15 

0.69 

19.45 

15.50 

15.35 


assuming Q > 1. The spectral- averaged quality factors of the conventional method (ICRP 26 in 
Anon. 1987; ICRP 60 in Anon. 1991) and proposed method (ICRU 40 in Anon. 1986) are shown 
in table 1 for the various isotopes produced in reactions. The proposed values are generally 
greater than the conventional values, except for the heavier fragments where the proposed 
values are substantially smaller. The conventional average quality factors show, generally, a 
weak isotope dependence, whereas the proposed average quality factors show a strong isotope 
dependence with neutron-rich isotopes being the most biologically damaging. 

Results and Discussion 

The dose equivalent per unit fluence of incident ions of charge Z and energy per nucleon 
E p are shown in figures 1 to 3. Figure 1 is based on current quality factors (ICRP 26 in 
Anon. 1987), figure 2 is based on newly proposed quality factors (ICRU 40 in Anon. 1986), and 
figure 3 is based on newly recommended values (ICRP 60 in Anon. 1991). The proton-induced 
fragmentation cross sections are taken from Wilson, Townsend, and Khan (1989). The proton 
cross sections are velocity-scaled according to the proposed factorization model of Lindstrom 
et al. (1975) as modified by Silberberg, Tsao, and Shapiro (1976). The limitations of this model, 
as discussed elsewhere (Wilson et al. 1984), do not concern us here because the 16 0 and 12 C ion 
beam data were used in the original derivation by Lindstrom et al. (1975), which was retained 
in subsequent modifications of Silberberg, Tsao, and Shapiro (1976), and adequately represent 
the 16 0 and 12 C data. The problems with this scaling model arise for nuclear fragmentation 
predictions far removed from projectile- target combinations used in fitting the model parameters. 
For example, there are no light- fragment data for iron fragmentation with which to fit the model 
(Wilson, Townsend, and Badavi 1987), but such experiments are currently in progress. The 
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Figure 1. Dose equivalent of various ion types includ- 
ing nuclear reactions for ICRP 2G quality factor 
(Anon. 1987). 


Figure 2. Dose equivalent of various ion types includ- 
ing nuclear reactions for ICRU 40 quality factor 
(Anon. 1986). 



Figure 3. Dose equivalent of various ion types including nuclear reactions for ICRP 60 quality factor (Anon. 1991). 


50-percent increase proposed for the 20 Ne data is within the uncertainty generally regarded for 
the Silberberg, Tsao, and Letaw (1983) parameterization. (See also Mathews (1983).) 

The nuclear contribution to the dose equivalent increases rapidly at the lowest energies as 
new channel thresholds are passed with increasing energy opening new reaction mechanisms. 
Only a small variation is seen between 20 and 300 MeV/nucleon and this is related to nuclear 
transparency (Townsend, Wilson, and Bidasaria 1982). New inelastic channels open above the 
pion production threshold and cause a rapid rise in dose equivalent above 300 MeV/nucleon. 
The fractional contribution of nuclear reaction effects is shown in table 2 for the two quality 
factors and ion types shown in figures 1 to 3 at 10 GeV/nucleon. The nuclear contribution to 
dose equivalent for C and heavier ions is less than 5 percent. Nuclear contributions for lighter 
ions can be substantial and as high as 70 percent. 

The average quality factors including nuclear reaction effects are shown in figures 4 to 6. The 
nuclear effects are clearly seen as the rise in average quality factor at high energies, especially 
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Table 2. Fractional Contribution of Nuclear Reactions to Total Dose 
Equivalent at 10 GeV/nucleon for ICRP 26, ICRU 40, and ICRP 60 Quality Factors. 



Fractional contribution for projectile- - 

Report 

'H 

4 He 

3 Li 

9 Be 

12 c 

28 Si 

56 Fe 

ICRP 26 

0.59 

0.43 

0.27 

0.17 

4.2xl0 -2 

2.5xl0 -3 

3.8x 1CT 4 

ICRU 40 

.70 

.51 

.34 

.21 

4.7 

2.1 

3.5 

ICRP 60 

.66 

.46 

.30 

.19 

4.6 

2.1 

3.5 



Figure 4. Average quality factor of various ion types 
including nuclear reactions for ICRP 26 quality 
factor (Anon. 1977). 



Figure 5. Average quality factor of various ion types 
including nuclear reactions for ICRU 40 quality 
factor (Anon. 1986). 



Figure 6. Average quality factor of various ion types including nuclear reactions for ICRP 60 quality factor 
(Anon. 1991). 


for the light ions. The increase in average quality factors for high-energy protons of ICRU 40 
(Anon. 1986) compared with the value obtained for ICRP 26 (Anon. 1987) is only 25 percent. 
Consequently, the earlier estimates of Alsmiller, Armstrong, and Coleman (1970), in which the 
quality factor of 20 is assumed for all nuclear fragments of mass greater than 1 amu (atomic 
mass unit), are expected to remain slightly conservative with respect to biological risk, even if 
the ICRU 40 (Anon. 1986) quality factor is enforced. 
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We now cite a few experiments that have measured relative biological effectiveness (RBE) 
for some biological end points. A genetic change in radiation-exposed seeds of maize produces a 
yellow-green streak in the leaves of the growing plant. The measured RBE for 28 GeV protons 
yielded an average of 4.4 (Smith 1967). The emergence rate of exposed brine shrimp was studied 
using 645 MeV and 9.2 GeV protons for which RBE values of 2.3 and 1.5 were reported (Gaubin 
et al. 1979). A third experiment on micronuclei induction in bean roots by 250 GeV protons 
found RBE’s of 1.8 to 2.1 (Diehl-Marshall and Bianchi 1981). A fourth experiment showed 
much more dramatic bone marrow depression in primates with 400 MeV protons compared with 
32 and 55 MeV protons of energy (Conklin and Hagan 1987). Although evidence exists of RBE 
being greater than unity because of target fragment effects, cancer induction in humans is not 
clearly represented by the quality factors presented herein. 

Effects on Cellular Track Models 

The cellular track model of Katz has been described extensively (Katz et al. 1971, 1972, 
1986). Here, we outline its basic concepts and consider the extension to the mixed-radiation 
field seen in space. The biological damage from passing ions is caused by delta-ray production. 
Cell damage is separated into a so-called grain-count regime (where damage occurs randomly 
along the ion path) and a track-width regime (where the damage is said to be distributed like a 
“hairy rope”). The response of the cells is described by four cellular parameters, two of which 
(m, the number of targets per cell, and D a , the characteristic X-ray dose) are extracted from 
the response of the cellular system to X-ray or 7 -ray irradiation. The other two parameters 
(£ 0 , interpreted as the cross-sectional area of the cell nucleus within which the damage sites 
are located, and k, a measure of size of the damage site) are found from survival measurements 
with track-segment irradiations by energetic-charged particles. The transition from the grain- 
count regime to the track-width regime takes place at Z* / Kp 2 4, where Z * is the effective 

charge and (3 is the velocity. The grain-count regime occurs at lower values of Z* and the 

track-width regime occurs at the higher values. 

To accommodate for the capacity of cells to accumulate sublethal damage, two modes of 
inactivation are identified: ion kill (intratrack) and gamma kill (intertrack). For cells damaged 
by the passage of a single ion, the ion-kill mode occurs. The fraction of cells damaged in the 
ion-kill mode is taken as P = £/E 0 , where £ is the single-particle inactivation cross section and 
P is the probability of the damage in the ion-kill mode. Cells not damaged in the ion-kill mode 
can be sublethally damaged by the delta rays from the passing ion and then inactivated in the 
gamma-kill mode by cumulative addition of sublethal damage caused by delta rays from other 
passing ions. The surviving fraction of a cellular population N (whose response parameters are 
m,D 0 ,E 0 , and k) after irradiation by a fluence of particles F (Katz et al. 1971) is written as 



N 

_ _ TTi X 7 r 7 

(15) 

where the ion-kill survivability is 

* = e - EF 

(16) 

and the gamma-kill survivability is 



7T 

= 1 - (l - e _Z V D °) m 

(17) 

The gamma-kill dose is 

b 

-j 

ii 

h-t 

1 

b 

b 

(18) 
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( 19 ) 


where D is the absorbed dose. The single-particle inactivation cross section is given by 


E = E 0 1 - exp (^-Z* 2 /kP 2 ^ j 
where the effective charge number is 

Z* = Z 1 — exp (-125/?/Z 2/3 )] 


( 20 ) 


In the track- width regime where P > 0.98, we take P = 1. 

For cell transformation, the fraction of transformed cells per surviving cell is 

T = 1 — (N'/N' 0 ) (21) 

where N' /N' 0 is the fraction of nontransformed cells, and a set of cellular response parameters 
for transformations m', D' 0 , E' 0 , and k! are used. The RBE at a given survival level is given by 

RBE = D x /D (22) 


w licit; ( . ■n 

D x = -D 0 In { 1 — [1 — (N/N 0 )] l/rn j (23) 

is the X-ray dose at which this level is obtained. Equations (15) through (23) represent the 
cellular track model for monoenergetic particles. Mixed-radiation fields have been considered 
previously in the Katz model. (See, for example, Katz, Sharma, and Homayoonfar 1972.) Next, 
we consider placing the model in terms of the particle fields described previously. 


Evaluation of the Katz Model 

The target fragmentation fields are found in closed form in terms of the collision density 
(Wilson 1977) because these ions are of relatively low energy. Away from any interfaces, the 
target fields are in a local equilibrium and may be written as 


(p a (x, Eq, ; Ej ) 


1 r°° do aj (E', Ej) 

S a (E a ) J Ea dE f 


(®> Dj ) 


dE' 


(24) 


where the subscript a labels the target fragment type, S a (E ) denotes the stopping power, and 
E a and Ej are in units of MeV. 

The particle fields of the projectiles and target fragments determine the level and type of 
radiological damage at the end point of interest. The relationship between the fields and the 
cellular response is now considered within the Katz cellular track model. 

The ion-kill term now contains a projectile term as well as a target fragment term as 


ro o 

(EF) = Ej (Ej) (pj (x, Ej) +J2 dEa 4 a Ea ; E j) 

a ^0 

whereas the corresponding gamma-kill dose becomes 

D~i — [l — Pj (Dj)] Sj {Pj) i x ’Pj) 

/•oo 

+ y dE a [1 - Pa (E a )j S Q ( E a ) <t> a (*, E a ; Ej) 

a J 0 


(25) 


(26) 
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Use of equations (24) and (25) allows one to define an effective cross section as 


E* (Ej) = Ej ( Ej ) + Y f°° iE 

a J 0 *(*{&<*) JE a 


,/ d^otj -Ej) 
dE / 


( 27 ) 


The first term of equation (27) is caused by the direct ionization of the media by the passing 
ion of type j. The second term results from target fragments produced in the media. 


Results and Discussion 

Katz (Waligorski, Sinclair, and Katz 1987) has obtained cellular parameters for survival and 
neoplastic transformations of C3H10T1/2 from the experiments of Yang et al. (1985) as given in 
table 3. We note the large uncertainties in the transformation data of Yang, which should lead to 
a similar uncertainty in the transformation parameters. Parameter sets were found from data for 
instantaneous and delayed plating of the cells after the irradiation. Here, only the delayed plating 
case is considered. General agreement with the measured RBE values (Waligorski, Sinclair, and 
Katz 1987) was found using these parameter sets. 


Table 3. Cellular Response Parameters for C3H10T1/2 Cells 



Values for response parameters- — 

Cell-damage type 

m 

D 0 , cGy 

E 0 , cm 2 

K 

Killing 

3 

280 

5.0 x 10 7 

750 

Transformation 

2 

26 000 

1.15 x 10~ 10 

750 


The single-particle inactivation cross sections neglecting target fragmentation of equation (27) 
are shown in figures 7 and 8 for cell death and transformation, respectively, as a function of the 
energy (in MeV/amu) of the passing ion. The target fragmentation contribution (the second 
term of eq. (27)) for protons has been evaluated as shown in figures 9 and 10. For protons, 
the effect of the target fragments (dashed line, the second term in eq. (27)) dominates over 
the proton direct ionization (dotted line) at high energy. For high-LET particles (low energy), 
the direct ionization dominates and target fragmentation effects become negligible. A simple 
scaling by s/A~j relates the proton target fragment term to ions of mass A y The resulting 
effective action cross sections for cell kill and transformation are plotted in figures 11 and 12, 
respectively. We note that the low-energy ®®Fe component of the GCR spectra extends into the 
track-width regime where E > E 0 , and it is not represented in the present calculation. The error 
introduced by the present calculation is small. 

The RBE is found in the Katz model through equation (22) and has a functional relationship 
to dose (Cucinotta et al. 1991a, 1991b; Katz and Cucinotta 1991) given in the low dose limit as 

RBE = °° (lif ) 1/m ° [(1/m) ~ ( 28 ) 

The implications of equation (28) on space exposure are discussed elsewhere (Cucinotta 
et al. 1991a). We might ask if the high-LET components of proton exposure show similar 
characteristics. The RBE for cell survival of Chinese hamster cells exposed to 160 MeV pro- 
tons was studied by Hall et al. (1978), and their results (dashed curve) are compared with the 
present model (solid curve) in figure 13. We note that the RBE, without accounting for target 
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Figure 7. Cell killing cross sections for several ions in Figure 8. Cell transformation cross sections for several 

C3H10T1/2 cells according to Katz model. Ioniza- ions in C3H10T1/2 cells according to Katz model, 

tion effects only. Ionization effects only. 




Proton energy, MeV Proton energy, MeV 

Figure 9. Cell killing cross sections including effects of Figure 10. Cell transformation cross sections including 

nuclear reaction effects for protons in C3H10T1/2 nuclear reaction effects for protons in C3H10T1/2 

cells according to Katz model. ce ^ s according to Katz model. 
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Figure 11. Cell killing cross sections including nuclear 
reaction effects for various ions in C3II10T1/2 cells. 



Figure 12. Cell transformation cross sections includ- 
ing nuclear reaction effects for various ions in 
C3H10T1/2 cells. 


fragmentation, remains near unity (dotted curve in fig. 13), and the rise in RBE observed at low 
doses is solely connected to the target fragments. (Compare with fig. 9.) 

Effects on Fluence- Related Risk Coefficients 


Harderian Gland Tumors 

The idea underlying the concept of risk per unit fluence has been used by researchers for ion 
exposure (Todd 1964; Curtis, Dye, and Sheldon 1965). The fluence-related risk coefficient Fj is 
defined as the probability of a given end point of interest (e.g., cancer) per unit fluence of type j 
particles passing through the organ (Curtis et ah 1992). A first estimate of Fj(Lj) can be found 
from the RBE values of Fry et al. (1985, 1986) as approximated (Wilson et al. 1991a, 1991b) by 


RBE - 0.95 + y{ 1 + 2e_L/14 ) l 1 ~ ex P (- a 2 L2 ~ «3 l3 )] ( 29 ) 

where a\ = 18 720, a .2 = 7.43 x 10 -G , and 03 = 1.14 x 10 -8 . (See fig. 14.) Using Curtis et al. 
(1992) and Wilson et al. (1991a, 1991b), the relation 


Fj ih) = 


BBE(Lj) Lj 
12.5D„ 


(30) 


represents the risk coefficient for direct ionization only because RBE was taken as unity for 60 Co 
gamma rays. 

In addition to the ionization caused directly by primary and high-energy secondary nuclei 
from fragmentation of the primary ions, the nuclei constituting biological tissue (i.e., the “target” 
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Proton RBE 



Figure 13. RBE predicted by Katz model for direct 
ionization and including target fragmentation in 
comparison with experimental results. 



L , keV/jim 

Figure 14. Relative biological effectiveness (RBE) for 
Harderian gland tumors as function of linear energy 
transfer. 


nuclei) will break up into lower energy and, in some cases, very highly ionizing target fragments. 
Target fragment fluences ), produced by a passing energetic ion of energy Ej } are given by 

* j ( e j) = Se ° j {El) fj ^ {Ei) dE 'i (31) 

where o’ J (E l ) is the macroscopic fragmentation cross section, fj(E Ej) is the energy distribution 
of the jth fragment, and 4>i{Ei) is the Auence of passing ions of energy Ej. The total prevalence 
at low exposure is given in terms of F % as 

coo 

p = Fi (Li) 4>i (Ei) + E 7 0 F i i L J (£;)] *3 {Ej) dEj (32) 


An effective F*(L) using equation (32) is deAned as 


F* (Li) = 


Fr 



dEj 


h \h fe )] 



(Ei) fj (. E',Ei ) 


(33) 


The effective risk coefficient F* is shown in figure 15 as a function of particle energy for 
representative charge components with the target fragment contributions shown separately. The 
nuclear contribution to the effective risk coefficient is reasonably approximated by 


2.0 x 10~ 3 


1.1 x 10~ 3 


p* ( T \ _ A QA J x x ^ I 

ni{ 1 ’ \ 1 + exp [- (E - 12) / 4.5] 1 + exp [- (E - 700) /363] J 


(34) 
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Risk cross section, pm 



Figure 15. Risk cross section as function of particle en- 
ergy for total contribution (ionization plus target 
fragmentation) and for target fragmentation contri- 
bution. 



Figure 16. Risk cross section F l (L l ) for prevalence of 
Harderian gland tumor as function of LET used in 
calculations. 



Figure 17. Relative biological effectiveness (RBE) for Harderian gland tumors including target fragment effects. 

where A is the ion mass, E is the ion energy in units of MeV/amu, and L, is the ion LET 
in keV/^m. The effective prevalence risk coefficient is shown in figure 16 for various ion 
components. The contribution from direct ionization is shown as the dashed curve. As yet, no 
known experimental data exist for proton exposures to test these nuclear fragment contributions 
to the prevalence risk coefficient. The effect of target fragmentation on RBE values is shown in 
figure 17 as the dashed extensions. 
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Cataract Index 


A similar formalism has been developed for the formation of stationary cataracts in young 
rabbits. The risk coefficient for a stationary cataract index from direct ionization (Shinn 
et al. 1991) is taken as 


F(L) = 


1.12 x 10~ 2 L + 4 x 10~ 4 L 2 
13.5 


( L < 170 keV//im) 
(L > 170 keV//zm) 


(35) 


and is shown in figure 18 in comparison with the data of Lett et al. (1988, 1989). Nuclear 
fragmentation contributions have been evaluated elsewhere (Shinn et al. 1991) and are shown in 
figure 19. The RBE of high-energy protons is calculated to be 1.4 in comparison with estimates 
from data for 2 GeV protons between 1 and 1.6 as estimated by Lett and Cox (from a private 
communication). 



L , keV/jjm Z., keV/gm 


Figure 18. Fluence-based risk coefficient for intermedi- Figure 19. Effect of target fragmentation on risk 

ate phase of cataract index. Data points arc from coefficient for intermediate cataract, 

left- to- right 60 Co 7 and 20 Ne, 40 Ar, and 56 Fe ions. 

Possible Effects in Human Cancer Response 

In an earlier section, we evaluated the effects of target fragments on values of dose equivalent 
for various ions and found important contributions for light ions. We now consider the more 
complicated task of estimating their effects on protracted exposure with energetic ions. The 
complications arise since human cancer risk coefficients are mainly known for low-LET-exposure 
data with single exposures. 

The lifetime excess cancer risk is given in NCRP 98 (Anon. 1989) for acute (2-day) exposure 
at age 35 as 

R a = 2.2 x 10 ~ 2 H [1 + {Hj 1.16)] (36) 

with a similar result for a 10-year protracted exposure beginning at age 35 as 

R p = 1.8 x 10 ~ 2 H (37) 

The dose equivalent is given as 

H = Q(L ) L [4>j (L) / 6.24] (38) 
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Q(L) 


where Q(L) is the quality factor, L is LET in keV//im, and is the type j particle 

fluence in number per fim 2 . Two quality factors are considered as recommended by ICRP 26 
(Anon. 1987) and ICRP 60 (Anon. 1991) as shown in figure 20. The corresponding fluence-based 
risk coefficients are 

F a {L,cj> J ) = Ra/ct>j(L) (39) 

and 

FpiL,^) = Rpl4>j[L) (40) 

and are shown in figure 21 for the ICRP 26 quality factor and in figure 22 for the ICRP 60 quality 
factor. The acute fluence risk coefficient is taken as the “initial slope” value of equation (36). 


According to the prescription of NCRP 98 
41 percent less than the acute exposure risk. 



Figure 20. Quality factors defined by ICRP 26 and 
ICRP 60 (Anon. 1991). 


non. 1989), the protracted exposure risk is 



Figure 21. Excess cancer risk according to ICRP 26 and 
NCRP 98 (Anon. 1989). 



Figure 22. Excess cancer risk according to ICRP 60 and NCRP 98 (Anon. 1989). 
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A correct understanding on the use of these risk coefficients is found in how data are used to 
derive dose- related risk coefficients and quality factors. The risk coefficients are dominated by 
acute low-LET exposures of exposed individuals in the nuclear weapons blasts of World War II. 
The dose rate reduction factor is derived mainly from controlled animal experiments using X-ray 
or 7-ray exposures. The quality factor is derived by concensus of learned individuals based on 
measured RBE’s for various biological systems and end points. The RBE is assumed to reach 
a maximum at low dose and/or dose rate on which quality factor is judged. One would assume 
that the quality factor is appropriate for low dose rate and that it generally overestimates for 
acute exposure. More recent evidence on life shortening in mice indicates a possible dose rate 
enhancement for high-LET exposures in the low dose region. This enhancement is indicated by 
the dashed curve in figures 21 and 22. Such enhancements are suggested to be the result of cell 
cycle effects (Rossi and Kellerer 1986) or, more simply, to be a property of the neutron response 
curve (Rossi 1981), or they may appear as repair-dependent phenomena. 



Figure 23. Excess cancer risk according to ICRP 60 
(Anon. 1991) and NCRP 98 (Anon. 1989) assuming 
a quality factor of 20 for fission neutrons. 


Fission neutrons 

y-rays 

Fractions 

24 



Figure 24. Life shortening in mice comparing exposure 
to neutrons with y-irradiation for 24 or 60 fractions 
and for a single exposure. (Data are taken from 
Thomson et al. 1981 1989.) 


Because life-shortening studies show that loss of life is enhanced by fractionation at high 
LET, we approximate the lifetime-cancer- risk cross section by twice the acute exposure values 
above 100 keV/pm. We assume that sparing is appropriate for the LET range for which Q = 1 
(i.e., L < 10 keV//im). A smooth connection is made between 10 and 100 keV//im from the 
protracted cross section at low LET to twice the acute cross section at high LET. This is shown 
as the dashed curves in figures 21 and 22. The corresponding response curves are shown in 
figure 23 for 7-ray and fission neutron exposure (Q « 20). The relative values are similar to 
the life-shortening data in mice (Thomson et al. 1981-1989) as shown for comparison in figure 24. 
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The fission neutron curve is the same as that expected for relativistic Ar ions. We have calculated 
the cross section for direct ionization (according to NCRP 98 (Anon. 1989)) given as 


( 3.5 x 10 _3 L 

F a = < 3.5 x 1(T 3 L (0.32L - 2.2) 
1 1.06 y/L 


(L < 10) } 

(10 < L < 100) V (41) 

(100 <L) ) 


with the corresponding protracted cross section assuming enhancement as 


{ 2.08 x 10 3 L 

2.08 x 10“ 3 L (0.75 + 0.025L) (0.32L - 2.2) 
2Ay/L 


(. L < 10) 
(10 < L < 100) 
(100 < L) 


(42) 


where L is LET in keV///m. Note that equation (42) corresponds to the dashed curve in figure 22. 
The target fragment contributions are reasonably approximated by 


K,a(E) = A 


0.4 


4.2 x 10 


-4 


9.6 x 10 


-4 


1 + exp [- (E - 12) / 4.5] + 1 + exp [- (E - 700) /363] } 


(43) 


and provide a reasonable approximation to the nuclear effects in acute exposure. 

In equation (43), A is the ion atomic weight and E is the kinetic energy in MeV/amu. The 
nuclear fragment effects on protracted exposure, assuming no dose rate enhancement, are also 
given by equation (43) if no repair occurs above 100 keV///m. Assuming dose rate enhancement, 
the nuclear target fragment contribution is given as 


, 0,1 f 8.4 x 10 -4 1.9 X] 

F n,p ( E )-A ■ | ! +exp f_( E _ 12 )/4.5] + 1 + exp [- (E 


700) /363 


(44) 


The fractional contribution of nuclear fragments (assuming dose rate enhancement) to the total- 
cancer-risk cross section is given in table 4. Significant target fragment contributions occur even 
for 100 MeV protons both to the acute exposure and especially to the protracted exposure 
assuming dose rate enhancement. The total cancer risk is dominated by target fragment 
contributions at high energy as shown in this table. The risk cross sections of various ions 
are shown in figures 25 and 26. 


Table 4. Fraction of Nuclear Fragmentation Contribution to 
Total-Cancer-Risk Cross Section for Protons 



100 MeV 

1000 MeV 

F* /F* 

* n,a/ x a 

0.14 

0.66 

F* /F* 

■* Tltf/ * p 

.35 

.85 


The excess-cancer-risk coefficients per unit fluence are used in connection with the galactic 
cosmic ray (GCR) transport code HZETRN (Wilson, Townsend, and Badavi 1987) to evaluate 
astronaut risk to 1 year of GCR exposure protracted over a 10-year career. The results are given 
in table 5 for various particles (neutrons, protons, alphas, L (for 3 < Z < 9), and H (for 10 < Z)), 
as well as for the total exposure. The values given are according to the ICRP 60 quality factor 
(Anon. 1991) using the BEIR V risk data. The risk coefficient labeled protracted in figure 22 
is used. The values of table 5 in parentheses assume dose rate enhancements at high LET and 
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10" 1 rr 



10 


-4 

10 ' 


10 u 


10 l 


L, keV/|im 



L, keV/pm 


Figure 25. Acute-exposure-risk cross section for pro- 
tons and alpha particles. 


Figure 26. Protracted-exposure-risk cross section for 
protons and alpha particles. 


Table 5. Excess Cancer Risk on Astronaut’s Career for 1-Year Total Exposure 
Behind Aluminum Shielding Thickness 

[Values in parentheses assume dose rate enhancement] 




Excess cancer risk, percent, for- 

- 


Aluminum shielding 







thickness, g/cm 2 

n 

P 

a 

L 

H 

Total 

0 

0 (0) 

0.27 (0.51) 

0.17 (0.24) 

0.28 (0.53) 

1.15 (2.90) 

1.81 (4.17) 

1 

0 ( .01) 

.21 ( .40) 

.07 ( .09) 

.23 ( .41) 

1.00 (2.54) 

1.52 (3.45) 

2 

.01 ( .02) 

.23 ( .42) 

.06 ( .09) 

.22 ( .39) 

.94 (2.34) 

1.45 (3.25) 

3 

.01 ( .03) 

.24 ( .43) 

.06 ( .09) 

.21 ( .36) 

.86 (2.17) 

1.38 (3.07) 

5 

.02 ( .05) 

.25 ( .44) 

.06 ( .09) 

.18 ( .32) 

.75 (1.87) 

1.26 (2.76) 

10 

.03 ( .08) 

.27 ( .45) 

.05 ( .08) 

.14 ( .24) 

.54 (I. 33 ) 

1.00 (2.18) 

15 

.04 ( .11) 

.28 ( .45) 

.05 ( .07) 

.11 ( .18) 

.40 ( .98) 

.87 (1.80) 

20 

.05 ( .13) 

.28 ( .45) 

.04 ( .06) 

.08 ( .14) 

.30 ( .74) 

.75 (1.52) 

30 

.06 ( .16) 

.27 ( .43) 

.03 ( .05) 

.05 ( .09) 

.18 ( .44) 

.59 (1.15) 

50 

.06 ( .18) 

.23 ( .36) 

.02 ( .03) 

.02 ( .04) 

.07 ( .17) 

.41 ( .77) 


correspond to the risk coefficient labeled protracted, enhanced in figure 22. Note that the shield 
effectiveness depends on the risk model used. An aluminum shield of 15 g/cm 2 reduces the 
protracted risk by a factor of 2, whereas the protracted, enhanced risk is reduced by a factor 
of 2.3. We expect a strong dependence on shield material selection since high-LET components 
are strongly dependent on material composition (Wilson et al. 1991). None of the risk coefficients 
discussed herein are recommended for design studies, but the present study may be helpful in 
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defining the import ance of the effects of dose rate enhancement in future recommended protection 
practices. 

Concluding Remarks 

In passing through a small region of tissue, the local energy deposit was shown to consist of 
direct ionization and excitation of electrons and energy given up to low-energy target fragments. 
The formalism was coupled to biological-related functions to evaluate the effects of target 
fragment contributions. In conventional risk assessment, the effective quality factor for high- 
energy protons is increased from its direct ionization value of 1 to 3.3 at 10 GeV. Experimental 
relative biological effectiveness (RBE) values between 1.4 and 4.4 have been observed for several 
biological systems for protons in the multi-GeV range. 

Target fragmentation effects on track-structure models reveal a significant increase in proton 
RBE at low r exposures. Such high RBE values result from the high efficiency of repair of X-ray- 
induced injury at low exposure and low dose rate. The predicted RBE is in good agreement 
with experiments for 160 MeV protons and Chinese hamster cells. Most of the injury from direct 
ionization of protons at low exposures is anticipated to be repaired at low exposures, and the 
main contribution to injury will result from target fragments for which repair is largely inhibited. 

Significant contributions of target fragmentation are predicted for cancer induction in the 
Harderian gland and lesser contributions to cataract formation. No repair-dependent models 
exist for these systems. A simple repair-inhibited model is given for human cancer induction 
within which the inhibited repair of target fragment contributions further shows their importance 
for protracted exposures. 

Fragmentation of target nuclei within biological systems has been shown to be an important 
contribution to light ion exposure. A primary concern is that direct ionization of light ions 
probably exhibits repair as an influence in radiation response, whereas the highly ionizing target 
fragments show little or no repair and may, in fact, exhibit enhanced risk at low exposure rates. 
Such factors have a potentially large impact on spacecraft shield design in future NASA missions. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
November 6, 1992 
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